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Abstract 

CO \ QCD sum rules are useful tools for studying the spectral properties of 

hadrons; however, assumptions underlying standard sum-rule analyses can 
lead to inconsistencies with known results of chiral perturbation theory. This 

■ possibility is demonstrated with QCD sum-rule extractions of the nucleon 
00 ■ mass and a term. In both cases, inconsistent chiral behavior leads to an 

uncertainty in the sum-rule predictions of ~ 100 MeV. 
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QCD sum rules are based on the analytic properties of current-current correlators and 
asymptotic freedom . Hadronic spectral properties are related via dispersion relations to 
the current- current correlator evaluated in terms of fundamental quark and gluon fields. This 
evaluation uses an operator product expansion (OPE), which, in turn, requires knowledge 
of various vacuum matrix elements — the condensates. In practice, applications of QCD 
sum rules require simple parameterizations of spectral functions and a truncated OPE. 
These simplifications necessarily lead to uncertainities in the extracted spectral parameters. 
Moreover, it is difficult to assess the magnitude of these uncertainties reliably. Accordingly, it 
is useful to test QCD sum-rule results against reliable physical and mathematical constraints. 
In this Letter, we will use well-known properties of hadronic spectra in the chiral limit to 
test the dependability of QCD sum-rule extractions of the nucleon mass and a term. 

It is generally believed that chiral perturbation theory accurately describes all low-energy 
observables of QCD in the limit of light current quark masses []3],|4]]. In the present context, 
however, we will only be interested in the leading nonanalytic behavior of observables as 
a function of m q , the average of the up and down current quark masses. It should be 
noted that this nonanalytic behavior can be determined without the full machinery of chiral 
perturbation theory — it depends only on the existence of dispersion relations and a pseudo- 
Goldstone pion. 

In principle, an exact evaluation of a current-current correlator in terms of QCD degrees 
of freedom must reproduce all of the hadronic physics, including the physics associated with 
the existence of a pseudo-Goldstone pion; therefore, an "exact" implementation of QCD 
sum rules must be consistent with chiral perturbation theory. As noted above, however, all 
practical implementations are not exact; thus practical applications of QCD sum rules are 
not necessarily consistent with chiral perturbation theory. By studying the magnitude of this 
inconsistency, one can estimate the scale of the uncertainities in QCD sum-rule extractions 
of spectral parameters. We find uncertainties of ~ 100 MeV for both the nucleon mass and 
a term. This implies that QCD sum rules, as generally used, are not suitable for evaluations 
of the a term. 

QCD sum-rule analyses of the nucleon mass || are based on the time-ordered correlation 
function Hn{q) defined by 



where |vac) is the physical nonperturbative vacuum state and t]n is an interpolating field 
with the quantum numbers of a nucleon. We take the interpolating field for the proton to 



where u a and d a are up and down quark fields (a is a color index), T denotes a transpose in 
Dirac space, and C is the charge conjugation matrix. The interpolating field for the neutron 
r) n is obtained by interchanging the up and down quark fields. Lorentz covariance and parity 
invariance imply that the Dirac structure of Hn{q) is of the form 0,^] 
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where II is the unit matrix in Dirac space, and 111 and U q are Lorentz scalar functions of q 2 
only. 

The key ingredients of QCD sum rules are analyticity and asymptotic freedom. The 
scalar functions ITj(s) (i = {l,q}) are analytic everywhere except on the positive s axis 
; thus an application of the Cauchy theorem yields 

[* dsWiWMUis) R =°° - I dsW t (s)Uf PE (s) , (4) 

JO J\s\=R 

where the weight functions Wi(s) are arbitrary entire functions; the contour integral on the 
right-hand side is evaluated in a counter-clockwise sense. We take R to be large, so we can 
evaluate the integral around the circle on the right-hand side of Eq. (|^) using the OPE M . 
The effects of a truncated OPE are considered below. The discontinuity An, is defined by 

AUi{s) = lim Ui(s + it) - U^s - it) . (5) 

Note that AIT is proportional to the spectral density. 

Although the OPE evaluation of Ilj(s) is only valid for large \s\ 0, it is useful to consider 
the analytic structure of rrp PE (s) for all s. One finds that np PE (s), like ITj(s), is analytic 
everywhere except on the positive s axis; therefore, we can introduce the following definition: 



r^^( S )Anp PE ( S ) R =°°- / dsW^U^is) . (6) 

JO J\s\=R 

Again, the contour integral on the right-hand side is evaluated in a counter-clockwise sense. 
By combining Eqs. (^) and (^|), we arrive at the following sum rule: 

/*oo /*oo 

/ dsW i (s)AU i (s)= / dsW^AUf^is) . (7) 
Jo Jo 

Note that the right-hand side of Eq. (|7|) is only defined in the context of Eq. (|]) . 

In order for the sum rules to be useful, one must choose the weighting function such 
that higher terms in the OPE are suppressed, since a truncated OPE is used in practical 
calculations. Such terms ultimately take the form of power corrections, which can be written 
schematically as 

n OPE (*) = E§- (8) 

Terms in the OPE proportional to higher- dimensional condensates correspond to higher 
values of n in Eq. (|]). A straightforward calculation then yields the following form for the 
OPE side of the sum rule: 



(9) 



One can use this result to test whether a particular weighting function efficiently suppresses 
higher terms in the OPE. In addition, one must choose the weighting function to suppress 
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the continuum contributions to the phenomenological side of the sum rule, which can be 
accomplished by choosing a W(s) that drops sufficiently rapidly with increasing s. Of course, 
there is no guarantee that both goals can be met simultaneously. 

The correlator contains contributions from the nucleon pole and from higher-mass (con- 
tinuum) states with the quantum numbers of a nucleon. Lorentz covariance implies that the 
correlator can be parameterized as |7],|| 

Hiv(g) = - x nj^ + u T(q) , (io) 

where specifies the strength of the coupling between the interpolating field and the 
physical nucleon state, and IT™ nt denotes the contributions from the continuum. Thus the 
scalar functions in the correlator are 

ni(s) = -A^^%r + nr(*), (n) 
n^) = -A^^ + nrw. (12) 

On the other hand, one can use the OPE to calculate the scalar functions in Eq. @. 
With the interpolating field in Eq. (0), one obtains PJlO| 



0)^M- S )«-9^fe)(W 2 > + - , (13) 
n ? PE («) = -6^1 M-*) " iabi ^XG^) 2 ) - |<gg> 2 + • • • , (14) 

where {qq) ~ —(225 MeV) 3 [qq = \(uu+dd)} and {(gG) 2 ) ~ 0.5 GeV 4 are the usual quark and 
gluon vacuum condensates . The ellipses denote contributions from higher-dimensional 
condensates, a s corrections, and contributions proportional to the current quark mass m q = 
\{m u + rrid). The latter contributions are not important for our purpose, since they yield 
terms that are analytic in m q ; we are mainly interested in terms that are nonanalytic in m q 
The usual QCD sum rules based on the Borel transform can be obtained by choosing 
Wi(s) = W q (s) = W(s) = e~ s / M , where M is known as the Borel mass fHH- Using the two 
Borel sum rules for the nucleon, one can readily solve for the nucleon mass: 



M, 



N 



/ dt wit) An° PE (t) - Am ont (t) 
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dsW(s) 



An? PE (s) - An c „ ont ( 



(15) 



It is here that one can see a potential limitation of the sum rules. 

We consider how the quantities in the sum rules depend on m q or, equivalently, m 2 = 

o o 

—2m q (qq) / ' j\. (All quantities of the form x denote the first term in the chiral expansion of 
x.) The leading terms in the chiral expansion of the nucleon mass are 

M N = M N + Arhl-^S^rhl + ... , (16) 

327T/2 
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where A is an unknown constant The expansion in Eq. ([16]) includes the leading non- 
analytic term in m q . In principle, the QCD sum-rule result for [Eq. (|15[) 1 reflects this 
expansion. 

In practice, however, it is difficult to obtain sum-rule results for the nucleon mass that 
are consistent with Eq. (|16j). According to Eq. (|1|), the 0(mfy contribution to the nucleon 
mass, which is nonanalytic in the current quark mass, must come from either the OPE or 
the continuum. We expect that it must be from the continuum; standard OPE calculations 
based on perturbation theory yield a sum of Wilson coefficients that are analytic in the 
current quark mass multiplied by vacuum condensates. We know of no vacuum condensate 
with an 0(m^) contribution. 

There are more serious difficulties. The leading chiral behavior of the quark condensate 

is 



11 



m 



2 



(qq) = (qq) 1 ~K ln-f + ■ ■ • , (17) 

where Mq is a mass parameter that is presumably ~ 1 GeV. We take Mq to be in the 
range 500-1500 MeV. The quark condensate given in Eq. fll7D is one in which perturbative 



contributions have been subtracted off fTTfl ; thus it is presumably suitable for use in the 
OPE. 

Thus AIlf PE contains an O(m^lnm^) chiral logarithm that must be eliminated some- 
where in the sum- rule expression for M^ [Eq. (|15|)], since it does not appear in the chiral 
perturbation theory expression for M N [Eq. ([!§)]■ It is clear that the contributions to the 
OPE correlator from higher- dimensional condensates will not eliminate this chiral logarithm, 
since they are multiplied by different powers of 1/q 2 ; therefore, the chiral logarithm must 
be canceled by contributions to the continuum. 

Thus we expect the continuum to play an important role in guaranteeing the correct 
chiral behavior of the nucleon mass as deduced by QCD sum rules. Standard treatments of 
QCD sum rules, however, assume a very simple form for the continuum 

Anr(3) = e( s -4)Anf PE ( s ), (is) 

where s l are continuum threshold parameters. We take = Sq = Sq. With this ansatz, 
Eq. (|T5j) becomes 

[ S0 dtW(t)ATl° PE (t) 

M N = ^0 • (19) 

/ dsW{s)AIl° PE {s) 
Jo H 

Such a simple treatment of the continuum will almost invariably lead to a spurious chiral 
logarithm in the resulting nucleon mass due to the 0(m\ In m 2 ) chiral logarithm in (qq). In 
addition, the needed 0{m^) contribution to the nucleon mass will almost certainly be lost. 
It is clear that one must describe the continuum contribution to the correlator quite exactly 
in order to obtain chiral behavior that is consistent with Eq. ([16]). 

What are the consequences of the spurious and missing nonanalytic terms in QCD sum- 
rule calculations of the nucleon mass? This question can be examined by an explicit Borel 
sum-rule calculation. Solving Eq.(|T9"D leads to the the following form: 
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M N = g(s ;M 2 )(qq) + -- - , 



(20) 



where the ellipsis denotes contributions from higher-dimensional condensates. We have 
shown that the leading term in ITp PE , which is proportional to {qq), leads to a spurious 
chiral logarithm contribution to the nucleon mass, given the usual simple ansatz for n™ nt . 
Contributions to the OPE correlator from higher-dimensional condensates do not resolve the 
problem, since they are multiplied by different powers of 1/q 2 ; therefore, for our purposes, 
we can neglect these contributions. Clearly Eq. fl20"D is not adequate for a precise determi- 
nation of the nucleon mass; however, it should be sufficient for determining the scale of the 
uncertainty in sum-rule estimates of due to spurious and missing nonanalytic terms. 

The nucleon mass, as determined by chiral perturbation theory, consists of unknown 
analytic (in m q ) and known nonanalytic (in m q ) contributions. In particular, from Eq. (|16|). 
one has 

M^ PT = analytic contributions — m\ + • • • . (21) 

It is believed that Eq. ( pT| ) gives the correct expansion of the nucleon mass [|]]. The QCD 
sum-rule estimate of the nucleon mass, as given in Eq. (|20|) , can also be written in terms of 
analytic and nonanalytic contributions: 

° ° 2 

M^cdsr = analytic contributions m tp> + ' ' ' > (22) 

32vr 2j2 M 



where we have used Eq. ( [17D and Eq. (|20|) in the chiral limit. While we recognize that the 
numerical value of the continuum threshold sq will vary with different choices of m q , we know 
of no realistic way to determine a chiral expansion of Sq analogous to the chiral expansions 
of M N and (qq) given in Eqs. (|T6| ) and (|T7|), respectively. It is difficult to describe the 



dependence of so, since it is an effective parameter in a simple continuum ansatz. Therefore, 
in deriving Eq. (^2j), we have neglected the m q dependence of so- We will return to this issue 
later. 

It is possible that numerical differences in the analytic contributions to Eqs. (|2l|) and 
(f2~2|) will partially compensate for the differences in the nonanalytic contributions to these 
equations; however, this would be purely fortuitous. We must assume that the differences in 
the nonanalytic terms lead to an intrinsic uncertainty SM^ in the QCD sum-rule estimate 
of the nucleon mass: 

5M N ~ -J^&Jln^l + -^ml ~ 65-110 MeV , (23) 
32tt»/J MS 32nfl 

where, in obtaining the numerical estimate, we have used physical values for all quantities. 
Thus the uncertainty is fairly modest if one is simply trying to determine a rough numerical 
value for the nucleon mass. 



Unfortunately, this is not the case with the nucleon a term, which can be defined as [12 



a N = 2m q J d 3 x ((N\qq\N) - (vac|gg|vac)) . (24) 
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The state vector \N) is for a nucleoli at rest. The most recent value of the a term is 
45 ± lOMeV @. 



By an application of the Hellmann-Feynman theorem, one can show [13 

o 

dM N 02 dM N 
aN = m ^ = m *lM[ ■ (25) 

The latter form is most useful for our purposes. Evaluating the derivative in Eq. ([25|) using 
a QCD sum-rule estimate of the nucleon mass leads to a QCD sum-rule estimate of the 



nucleon a term. This approach is equivalent to that used in Ref. fl4|| , which is based on 
evaluating the nucleon correlator in the presence of a background scalar field. 
In analogy to Eq. ([H]) , we can write the chiral expansion of the a term as 

a N = analytic contributions — m n + • • • ■ (26) 

Likewise, from Eq. (|22|), one obtains the following expansion of the QCD sum-rule estimate 
of the a term: 



3M N no , ml 



^.qcdsr _ ana jytj c contributions %- mi In — ^ + • • • . (27) 

32vr 2 /2 Mq 

Thus QCD sum-rule estimates of the nucleon a term that are based on a simple ansatz for 
the continuum have an uncertainty 5o~n given by 

5a N ~ -^^ml In + ~ 70-115 MeV . (28) 

32n 2 fl M o 64vr/2 

Unlike the case of the nucleon mass, the QCD sum-rule determination of the nucleon a term 
is completely unsatisfactory. The uncertainty in the QCD sum-rule estimate of the a term 
associated with chiral physics is about two times the value of the a term itself. 

At this stage, it is worth returning to the issue of the m q dependence of s . One might 
hope that a nontrivial dependence exists and that such a dependence can eliminate or reduce 
the mismatch between chiral perturbation theory and the QCD sum-rule description. For a 
given choice of weighting function, one can, in fact, fabricate an m q dependence of so that 
does eliminate the discrepancy. However, such a solution is completely arbitrary; we have 
no fundamental reason to believe that this occurs. Moreover, this solution to the problem 
is unsatisfactory, since the m q dependence that eliminates the mismatch depends on the 
choice of weighting function. Furthermore, in actual implementations of QCD sum rules, 
So is determined via some optimization procedure. Such a procedure is extremely unlikely 
to give rise to the m q dependence needed to achieve consistency with the correct chiral 
behavior. 

In conclusion, we have shown intrinsic uncertainities of ~ 100 MeV in QCD sum-rule 
extractions of the nucleon mass and a term associated with inconsistencies between the 
sum rules and known chiral behavior. As a practical matter, this suggests that attempts 
to improve QCD sum-rule calculations for the nucleon via an improved treatment of the 
OPE or a more judicious choice of weighting function are of limited value in reducing the 
uncertainties. We see from the present study that a 100-MeV uncertainty will persist until 
the inadequacies of the continuum ansatz are addressed. 
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